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ABSTRACT 

We present a panoramic narrow-band imaging survey of [O II] emitters in and around 
the C1G J0218. 3-0510 cluster at z=1.62 with Suprime-Cam on Subaru telescope. 352 
[O n] emitters were identified on the basis of narrow-band excesses and photometric 
rcdshifts. We discovered a huge filamentary structure with some clumps traced by 
[O n] emitters and found that the C1G J0218. 3-0510 cluster is embedded in an even 
larger super-structure than the one reported previously. 31 [O n] emitters were spectro- 
scopically confirmed with the detection of Ha and/or [O in] emission lines by FMOS 
observations. In the high density regions such as cluster core and clumps, star-forming 
[O n] emitters show a high overdensity by a factor of more than 10 compared to the 
field region. Interestingly, the relative fraction of [O n] emitters in photo- z selected 
sample does not depend significantly on the local density. Although the star forma- 
tion activity is very high even in the cluster core, some massive quiescent galaxies also 
exits at the same time. Furthermore, the properties of the individual [O n] emitters, 
such as star formation rates, stellar masses and specific star formation rates, do not 
show a significant dependence on the local density, either. Such lack of e nvironmental 
dependence is consistent with our earlier result by lHavashi et al.l (1201 it ) on a z = 1.5 
cluster and its surrounding region. The fact that the star-forming activity of galaxies 
in the cluster core is as high as that in the field at z ~ 1.6 may suggest that the 
star-forming galaxies are probably just in a transition phase from a starburst mode 
to a quiescent mode, and are thus showing comparable level of star formation rates to 
those in lower density environments. We may be witnessing the start of the reversal 
of the local SFR-dcnsity relation due to the "biased" galaxy formation and evolution 
in high density regions at high this redshift, beyond which massive galaxies would be 
forming vigorously in a more biased way in proto-cluster cores. 

Key words: large-scale structure of Universe - galaxies: clusters: individual: C1G 
J0218. 3-0510 - galaxies: evolution 
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It is widely known that the formation and evolution of 
galaxies strongly depend on their surrounding environ- 
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ments. Galaxy clusters are most massive object and the 
densest structures in the universe. As such, clusters and 
their surrounding regions serve as ideal sites for study- 
ing the roles of galaxy environment on galaxy formation 
and evolution. In the local universe, it is known that 
the star formation has already ceased long time ago in 
high density regions such as clusters and most of the 
galaxies therein are only passively evolving since then, 
while galaxies in the general field are still actively form- 



ing stars even at the present day (e.g ., Lewis I. et alj |2002| : 



I Gomez et al.l 120031 ; iKauffmann et aj f l2004h . Such environ- 
mental variation has been well recognized as a sharp cor- 
relation between gal axy morphology and local number den- 
sity of galaxies (e.g., Dres sier 198(1 Postman fc Gellerll 19841 ; 



IWhitmore et al] 1 19931: iDressler et alj Il997t ). In clusters at 



Muzzin et al] (|201lh found that the specific star for- 



mation rate (SSFR) of star-forming galaxies is independent 
of environment at fixed stellar mass but the fraction of star- 
forming galaxies is decreased in high density environment. 
This suggests that the environmental-quenching timescale 
would be rapid. Recent observations are discovering very 
distant clusters at z ~ 1.5 — 2.0 dBrodwin et al. 20 111: 



Fassbender et all 2011 ; IPapovich et alT 2010l ; Tanaka et al] 



20ld : iHenrv et al] l20ld : iGobat et alj l201ll ). It is crucial 



to observe clusters at such frontier redshifts so as to di- 
rectly witness and understand the early evolution of galax- 
ies and the physical processes which are responsible for the 
strong environmental dependence of galaxy properties at 
later times. 

To survey star-forming galaxies with a narrow-band fil- 
ter, that captures nebular emission lines such as Ha and 
[O 11] , is very efficient and effective for studying cluster 
galaxies because we can construct a large sample of star- 
forming galaxies without any expensive spectroscopic obser- 
vation. The correlation between the star formation activity, 
traced by a nebular emission, and the environment has al- 
ready been investigated by many authors rathe r intensively. 
For fi eld environment and groups at z =0.84, ISobral et al] 
(|201ll ) have constructed a large sample of Ha emitters with a 
narrow-band filter as a part of the High-z Emission Line Sur- 
vey (HiZELS), and found that the fraction of star- forming 
galaxies falls sharply as a function of local number density of 
photo- z sample from about 40 percent in the field to almost 
zero in rich groups. For clusters at z ~ 0.5-1.5, our previ- 
ous studies on Ha emitters in CL0939+4713 (Abell 851) at 
z = 0.41 and in RX J1716. 4+6708 at z = 0.81 show that the 
fraction of emitters decreases towards cluster central region, 
and the star formation activity peaks in the intermediate 
density re gions such as groups and fi laments at the cluster 
outskirts (|Kovama et alj l201ll . l20ld ). On the other hand, 
the studies in higher redshift clusters (z ~ 1.5) have pro- 
duced two different results. While there are no Ha emitters 
within a radius of 200 kpc from the centre of t he most mas- 
sive c luster XMMU J2235. 3-2557 at z = 1.39 (|Bauer et al] 
I2011D, the other clus ter XMMXCS 2215.9-1738 at z = 1.46 
( Havashi etHfllbOld ) shows much higher star formation ac- 
tivity as traced by [O 11] emitters even in the cluster core. 
One possible interpretation would be that, in a more massive 
system galaxy evolution is accelerated and star formation is 
truncated at an earlier epoch. However, there is an caveat 
that some of the [O 11] emitters in the cluster core may be 
contaminated by AGN. 



To obtain more general picture on the star formation 
history in clusters, we need to study many clusters at each 
epoch. In fact, we are currently conducting the MAHALO- 
Subaru project (MApping HAlpha and Lines of Oxygen with 
Subaru: Kodama et al. in prep), which is a large and sys- 
tematic narrow-band survey of star forming galaxies in many 
clusters at z ~ 1.5-3.0 and in the field. This paper reports 
one of the initial results of the project. This survey will pro- 
vide us with the dependence of star formation history on 
environment. 

We comment here, however, that if we only target clus- 
ters that have already formed and matured at each epoch, 
the results may be biased towards apparently weaker evolu- 
tion. This difficult problem should be tackled by targeting 
lower-mass systems as well as the rich clusters. By tracing 
the surrounding large scale structures around rich clusters, 
however, we can automatically include many smaller mass 
groups and /or filaments embedded in the structures. This 
is another good advantage of panoramic studies of distant 
clusters like ours. 

In this paper, we report a narrow-band survey of 
[O 11] emitters in and around the spectroscopically con- 
firmed dis tant cluster C1G J0218.3-0510 (IRC 021 8- A) at 
2 = 1.62 jPapovich et alj|20ld ; [Tanaka et al] 120101 ) and is 
in Subaru/XMM - Newton Deep Survey Field (SXDS; 
iFurusawa et al] 120081 ). The coordinates of this cluster is 
a = 2 h 18 m 21.3 s ,5 = -5°10'27"(J2000), derived from the 
centroid of the sources selected by Spitzer IRAC colour 
in this overdensity. We assume cosmological parameters of 
H = 70 km s _1 Mpc -1 , Sl M = 0.3, and fl A = 0.7, and 
adopt AB magnitudes throughout this paper. At z = 1.62, 
l"corresponds to 8.47 kpc in a physical distance and to 22.20 
kpc in a comoving distance. 



2 OBSERVATIONS AND DATA 
2.1 Imaging data 

To identify line emitters, we used a narrow-band filter 
(NB973) and a medium-band filter (zr). The NB973 filter 
(Ac = 9755A, AA FWHM = 202A) can detect an [O 11] line 
emission at 1.590 < z < 1.644 and z R filter (A c = 9860 A, 
AAfwhm = 590A) can measure a continuum flux at same 
wavelength range as the NB973 filter. These filters enable 
us to sample most of the star-forming galaxies associated 
to the cluster at z = 1.62 as shown in Figure [1] The imag- 
ing data with NB973 filter wa s obtained i n an area between 
the SXDS-C and SXDS-S by lOta et ail (|20ld ). The seeing 
size in the combined image is 1.0" and the 5<r limiting mag- 
nitude is 25.4 in a 2.0" diameter aperture. The details of 
the observations and the data reduction are described in 
lOta et all (|20ld ). For the zr data, we conducted a new ob- 
servation with Subar u Prime Focus Camera (Suprime-Cam; 
iMivazaki et alj|2002T l on the Subaru telescope in October 
2010. The weather was good during our observation, and 
the sky conditions were photometric with the seeing of 0.5 
- 0.7". Data reduc tion was carried out in a normal man- 
ner with SDFRED (|Yagi et al.ll2002i ; IOuchi et alj|2004 ). The 
frames for the total of 327 minutes are integrated in the fi- 
nal image. The field coverage is 830 arcmin 2 , corresponds 
to 1.4 x 10 5 Mpc 3 in the survey volume. The 5a limiting 
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Figure 1. A transmission curve (solid line) of NB973 filter and 
the redshift distribution (black histogram) of spectroscopically 
confir med galaxies, including the cluster memb ers in the litera- 
tures l|Papovich et al]|2010l ; iTanaka et al.ll2010h . The [O II] emit- 
ters, whose redshifts are derived from our FMOS observation, are 
shown by red histogram. 

magnitude is 25.3 in a 2.0" diameter aperture. Point spread 
function (PSF) of the zr image was degraded to 1.0" to 
match the quality of the NB973 image. 

2.2 Photometric catalogue 

In the SXDS field, many multi-wavelength data are al- 
ready available. In this work, we used a total of 12 multi- 
band data including the NB973 and ZR-h&nd ones. Deep 
optical broadband images (B ,V, R,i' , z ') were taken by 
the Suprime-Cam |Furusawa et al.ll200ct ). For near- infrared 
wave length, the UKIDSS Ultra Deep Survey (UKIDSS- 
uds|l awrence et al.ll2007f ) has been conducted in this re- 
gion, and we used the DR8 data (J,H,K). T he photomet- 
ric sy stem and calibration are d escribed by iHewett et al.l 
l|200fj ) and iHodgkin etafl (|2009l ), respectively. From 10- 
band (NB973,B, V, R, i' , z' , zr, J, H, K) i mages, we made a 
phot ometric catalogue using SExtractor (| Bert in fc Arnouta 
1 19961 ). PSFs of all the images were matched to 1.0"and the 
positions were matched with respect to the NB973 image. A 
source detection was performed on the NB973 image. The 
extraction criteria were at least 9 pixels with fluxes above 
2a level, where la is the sky noise of the image. The sources 
fainter than 5a limiting magnitude in NB973 were rejected. 
Photometries in all the images were carried out using the 
double image mode. An aperture magnitude within a di- 
ameter of 2"was used to derive a colour and Kron magn i- 
tude was used to calculate a total magnitude (|Kronlll980D . 
For mid-infrared wavelength, we used IRAC catalogue (only 
3.6 fim and 4.5 fim) for a Spitzer Public legacy survey of 
the UKIDSS-UDS (SpUDS; PI: J. Dunlop). From our cata- 
logue, we identified sources that have an association in the 
IRAC catalogue within a 1.0 arcsec radius. Moreover, not 
only the NB973-detected catalogue but also the K-detected 
one was created above 5a in the same manner when we select 



passively evolving galaxies later in section 3.3. The NB973- 
detected and the K-detected catalogues include 39229 ob- 
jects and 40786 objects, respectively, of which 33242 objects 
are in common. 

2.3 Spectroscopy 

Near-infrared spectroscopy of our photometrically identified 
[O n] emitter candidates wa s conducted with Fibe r Multi 
Object Spectrograph (FMOS; lKi"mura M. et~al]|2010l ) on the 
Subaru telescope in January 2012. Observations were made 
with low-resolution (LR) mode in IRS1 and high resolu- 
tion (HR) mode at J-long and H-long bands in IRS2. Low- 
resolution, J-long and H-long modes can cover the wave- 
length ranges of 1.0-1.8 /xm, 1.12-1.35 /xm and 1.59-1.80 
/j,m, respectively. Cross-beam switching mode was adopted 
and the total integration times were 375 minutes in LR, 
60 minutes in J-long and 300 minutes in the H-long mode. 
The seeing was typically 1" in R-band during observations. 
Red uction was carried out with the FMOS pipeline FIBRE- 
pac |lwamuro et al.l 120111 ) . The flux calibration was done 
using some F,G or K-type stars from the Two Micron All 
Sky Survey (2MASS; ISkrutskie et al.|[2006l ). observed simul- 
taneously. The flux loss was estimated to roughly 50% for 
a point source and this might be caused by various factors 
(e.g. weather condition, instrument focus, and position error 
in the catalogue). In this paper, only spectroscopic redshifts, 
derived from fitting emission lines, are used in order to esti- 
mate the accuracy of the photometric redshift, and the level 
of completeness and contamination. Other analyses such as 
line fluxes, line widths and the line ratios will be presented 
in detail in our forthcoming paper. A lso, we used th e cata - 
log ues of spectroscopic redshifts from lSimpson et al.l (|2012l ) 
and lSmail et al.1 (|2008l ). 



3 TARGET SELECTION 

In order to investigate galaxy evolution, we need to sample 
both actively star-forming galaxies and passively evolving 
galaxies at the same epoch. As noted in section [T] the for- 
mer can be selected by the presence of nebular emission 
lines such as [O n] . The latter population can be selected 
by utilizing the photometric redshift and the colour-colour 
diagram. Using the multi-band data in section 2, a total of 
352 [O n] emitters at z ~ 1.62 have been identified over a 
830 arcmin 2 area. Also, we have constructed a photo- z se- 
lected sample with M Bta r > 10 10 M at 1.56 < z phot < 1.68, 
which includes 132 [O n] emitters and 259 quiescent galaxies 
at z ~ 1.6. 

3.1 [O n] emitters 

We can identify emission line galaxies on the basis of the ex- 
cesses of NB973 fluxes over zr fluxes. NB973-detected cata- 
logue created in section 2.2 is used. Figure [2] shows the zr — 
NB973 colour-magnitude diagram for all the objects satisfy- 
ing NB973 > 5a. Note that we made a correction of +0.11 in 
the zr — NB973 colour because the effective wavelengths of 
two filters are not exactly same. For faint objects (< 2a) in 
z/j-band, 2a value was used to calculate colour. If an strong 
emission line comes into the zr band but not on the NB973 
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Figure 2. z^ - NB973 colour-magnitude diagram to select 
[O II] emitter candidates at z=1.62. Gray dots show all objects 
with NB973 > 5<r. Black and red filled circles are NB973 emitters 
and [O II] emitters which satisfy the criteria in section 3.1. 



filter (e.g., at z = 1.58 or z = 1.66), such object may be de- 
tected as false-absorbers. Also, because the effective wave- 
length of NB973 filter is slightly shorter than that of the 
zr filte r, very red objects m ay have negative colours in zr- 
NB973. [Bunker et al.l (|l995T l have defined the significance of 
the excess in the narrow-band by the parameter £ with tak- 
ing account of the fact that the fainter the narrow-band flux 
is, the larger the photometric error is. S =2.5 was adopted a s 
our first criteria to select NB emitters l|Sobral et al.|[2009T ). 
Since the standard deviation of zr — NB973 is 0.03 in the 
range of 18<NB973<21, the colour excess of z R - NB973 
> 0.1, that is corresponds to 3cr for the bright objects, was 
also adopted. These criteria correspond to the limiting line 
flux of 1.8 xlO" 17 ergs" 1 cm" 2 , that is 4-5 M Q yr _1 in dust- 
u ncorrected star fo rmation rate (SFR) using the calibration 
of lKennicuttl (|l998h . and to the equivalent width of about 30 
A in the observed frame. Based on first criteria, 969 objects 
are identified as NB973 emitters. 

The NB973 filter can pick out not only [O n] lines 
at z = 1.62 but some other lines at different redshifts: 
Ly a at z = 7.02, H/3 at z = 1.01, [O in] at z = 
0.95 and Ha line at z = 0.49. To discriminate between 
[O n] emitters and those other line emitters, we utilize 
photometric redshifts for those objects that have satis- 
fied our criteria. Ph otometric redshifts w ere determined 
with the EAZY code (|Brammer et al.ll2008fi . 11 band data 
(B, V, R, i' , z' , zr, J, H, K, 3.6/xm, 4.5^im) are used for SED 
fitting. Capturing Balmer/4000A break, z/j-band should im- 
prove the accuracy of photometric redshift for galaxies at 
z ~1. 4-2.0. Because z.R-band can be substantially affected 
by the presence of our target emission line itself, the emis- 
sion line flux and the emission line subtracted continuum 
flux density are estimated from the Zij-band and NB973- 
band fluxes as follows; 



Figure 3. The distribution of photometric redshift for NB973 
emitters. Red dash lines indicate the expected redshift of emission 
lines that can fall onto the NB973 filter (Ha at z ~ 0.49, [O III] at 
z ~ 0.95, H/3 at z ~ 1.01 and [O Ii] at z ~ 1.62). Blue lines 
indicate the photometric redshift ranges (1.4 < z p i 10 t < 1-9) to 
define [O Ii] emitters. 
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where 7<ii ne is the emission line flux, / denotes a flux den- 
sity, and A indicates a FWHM of a filter. The emission line 
subtracted flux density in z,R-band is used for SED fitting. 
Figure [3] is the distribution of the photometric redshift for 
NB973 emitters. We can clearly recognize three peaks in the 
photometric redshift distribution (i.e., z ~ 0.5, 1, and 1.6). 
These redshift peaks neatly correspond to Ha, [O ill] or H/3, 
and [O ii] respectively, assuring the detection of line emit- 
ters as expected. This also indicates that our first criterion 
of £=2.5 and zr — NB973 > 0.1 is effective in identifying se- 
cure line emitters. Out of NB973 emitters, we down selected 
352 [O ii] emitters that fall between 1.4 < z p hot < 1.9, as sec- 
ond criterion, to distinguish between our target [O n] emit- 
ters at z = 1.62 and other unwanted line emitters at other 
redshifts. 

While most of the [O n] emitters are likely to be star- 
forming galaxies, there is a possibility that some of the 
[O ii] emission lines are originated from central active galac- 
tic nuclei (AGN) rather than star forming regions. In par- 
ticular, some authors reported that "red" [O n] emitters 
tend to be contamin a ted by AGN activitie s jYan et al.ll2006l; 
Lemaux etal] |2010t lHavashi et"aH l201ll : iTanaka M. et all 
201lT ). Alternatively, red emitters are dusty starburst galax- 
ies which are reddened by a large amount of dust extinc- 
tion. In any case, these red emitters are interesting objects 
which may be key populations to understanding galaxy evo- 
lution in dense environment. However, we can not disentan- 
gle these two possibilities at this stage without spectroscopy. 
With spectroscopy, we should be able to distinguish between 
AGNs and dusty starbursts by measuring emission line ra- 
tios of [O ml/Hfi and [N h]/Hq and plot them on th e BPT 
diagram (|Baldwin et al.lll98ll ; lkauffmann et aLlbOOSt ). Here, 
we defined red emitters with the z — J versus J — K colour- 
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colour diagram (hereafter zJK diagram) and treated only 
the [O n] emitters that are outside the quiescent zone de- 
fined in section 3.3 as star-forming galaxies. Also, our cri- 
teria of [O n] em itters select three X-ray detected objects 
|Ueda et alj|200sh . of which two are located at different red- 
shifts. Therefore, the remaining one is not used in our anal- 
ysis due to some uncertainty. Eventually, we have identified 
totally 340 star-forming [O n] emitters and 12 red emitters 
associated to the cluster at z — 1.62. 

3.2 Spectroscopic follow-up for [O n] emitters 

At z = 1.62, [O in] and Ha emission lines shift to 1.31 fira 
and 1.72 fj,m, respectively. From the NB973 emitter sample 
with F([0 n])> 3.0 x 10~ 17 erg s^ 1 cm ~ 2 identified in the 
previous section, 46 (IRS1, LR mode) and 39 objects (IRS2, 
HR mode) were observed with FMOS. Ha emission lines 
are detected for 4 objects in LR mode and for 24 in H-long. 
[O in] emission lines are detected for 7 objects in LR mode 
and for 4 in J-long. The sensitivity of the LR mode is not 
good enough to detect emission lines due to the large flux 
loss. If only the H-long mode is considered, the detection 
rate of emission lines is about 60 percent in total and 80 
percent for bright [O n] emitters with F([0 n])> 7.0 x 10~ 17 
erg s _1 cm ~ 2 . These value are really good since 20-25 per- 
cent of the wavelength coverage is masked in the FMOS 
spectrograph due to OH-airglow lines. Figure U shows some 
examples of Ha spectra of our [O n] emitters. Spectroscopic 
redshifts are derived by gaussian fitting with free param- 
eters of redshift, line width, and flux densities of Ha and 
[N n]A6584 or [O m]A4959, A5007. As a result, 31 spectro- 
scopic redshifts are obtained in total from the FMOS data, 
since there are 8 objects whose Ha and [O ill] lines are both 
detected. There are a total of 42 spectroscopically confirmed 
galaxies with emission lines, including 11 spectroscopic red- 
shifts with a emission line in the literatures. 40 out of them 
were selected by the criteria in the previous section. On the 
other hand, there were no objects with a different redshift 
in our [O n] emitter sample. 

3.3 Quiescent galaxies 

We also need to select passively evolving galaxies located 
at 2=1.62 to quantify the averaged star formation activ- 
ity. First, photometric redshifts were measured for all the 
objects in the K-band detected catalogue created in sec- 
tion 2. X-ray detected objects were rejected because the 
accuracy of photometric redshift is so bad that could in- 
crease contaminations. In FigureO photometric redshifts are 
plotted against spectroscopic redshifts for confirmed galax- 
ies. The standard deviation of (z S p CC — 2phot) is 0.05 at 
1.590 < 2 sp cc < 1.644. Our photometric redshifts are fairly 
good for them. We selected only the objects that fall within 
a narrow redshift interval of 1.56 < z p hot < 1-68 since the 
K-band selected galaxies can include contaminations at con- 
tiguous redshift unlike the NB973 emitters. Such a stringent 
criterion of photometric redshift range is reasonable because 
our photometric redshifts are good enough to recover most 
of the spectroscopically confirmed cluster members. In fact 
the photometric redshift distribution of the [O n] emitters 
shows a narrow concentration at z ~ 1.6 as shown in Figure 

El 
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Figure 5. The spectroscopic redshift versus the photometric red- 
shift for galaxies with 0.8 < z S pcc < 1.8. Red line is z p ^ ot = z Bpcc . 
Blue dashed lines show the criterion of the photo- z selected sam- 
ple (1.56 < 2 p hot < 1.68). 




J-K 

Figure 6. The z — J versus J — K colour-colour (zJK) dia- 
gram of the photo-z selected galaxies (black dots). Blue circles 
show the [O II] emitter samples. Red dashed and green dashed 
lines shows a model track of a passive gal axy and a star-forming 
galaxy, respectively jKodama et alJll99Sf) . The squares indicate 
the model colour at z = 0, z = 1.6 and z = 2.5. The arrow show s 
the dust reddening vector of E(B-V)=0.2 jCalzetti et al.ll200tj) . 
The solid lines show the adopted criteria to separate between 
quiescent galaxies and star-forming galaxies at z = 1.6. 

We have constructed a mass-limited sample with the 
stellar mass M sta r > 10 10 M© on the basis of photometric 
redshifts. Stellar masses are estimated from total K-band 
magnitudes, K total , and z — K colours, using the ratio be- 
tween mass-to-luminosity ratio in K-band and z-K colour 
of the population sy nthesis bulge-disk composite models 
l|Kodama et al.|[l999l ) after it is scaled to the Salpeter IMF 
for consistency. At z = 1.6, we used the following equations: 
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Figure 4. Examples of Hq spectra for [O II] emitters. Red line and blue dotted line shows a observed spectra and a gaussian fitting 
spectra, respectively. Vertical dashed lines indicate the redshifted wavelength of Ha and [N ll]A6584. Gray zones indicate the positions 
of the OH-mask. 
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MAfstar/lO 11 ) = -0.4(Jf* 



21.64) 



AlogM s ta r = 0.1 - 1.12exp[-0.82(z' - K)] 



(3) 
(4) 



The amount of contamination and that of completeness of 
the sample are estimated based on the spectroscopically con- 
firmed galaxies. There are 32 galaxies with M s tar > 10 10 Mq 
which are located within 1.590 < z spcc < 1.644. Since 6 spec- 
troscopically confirmed cluster member galaxies do not fall 
within the photometric redshift range of 1.56 < z p hot < 1.68, 
the completeness can be estimated to be 81%. We include 
them in our sample. Also, because our sample include three 
objects at different spectroscopic redshift, from the clus- 
ter redshift, the contamination can be estimated to be 
10%. However, in order to accurately estimate them, we 
need to make more spectroscopic follow-up observations. 132 
[O n] emitters with M sta r > 10 10 Mo satisfy our photomet- 
ric redshift selection (1.56 < Zphot < 1.68). After excluding 
these [O n] emitters from our photo-z selected sample, we 
constructed a non-emitter sample of 1174 galaxies. 

It is known that cluster galaxies exhibit a conspicuous 
"red sequence" on the colour-magnitude diagram which are 
mainly composed of passively evolving galaxies. The red se- 
quence is well recogni zed even in distant clusters out to z ~ 1 
l|Kodama et aljfl998t ) and it persists at the brigh t end even 
in proto-clusters at z ~ 2 |Kodama et al.l [20071 ). The red 
sequence consists of not only passive galaxies but also dusty 
star-forming galaxies, and we use a zJK diagram to separate 
them in the same manner as in I Williams et "all (|2009l ). Figure 
[6] shows the zJK diagram for th e photo- z select e d gala xies. 
Using the model colour tra cks of Kodama et alj (|l998T ) and 
the dust reddening vector (|Calzetti et al.l 2000] ). we defined 
the quiescent galaxies at z = 1.6 using the three criteria: 
2 - J > 1.3, J-K < 1.6 and z - J > 0.64( J - K) + 0.8. In 
total, 259 quiescent galaxies have been selected from photo- 
z selected sample. On the other hand, there are a number 
of blue galaxies (915) in photo- z selected sample. It is likely 
that either their [O n] emissions are too weak to be detected 
as the narrow-band excess or they are outside of the redshift 
range of the NB973 filter (z < 1.590 or z > 1.644). Since the 
contamination of the latter objects are probably large, we 
do not use those blue galaxies in our analyses. Three qui- 
escent galaxies have been spectroscopically confirmed and a 
spectral fit of the brightest member yields an age of 1.8 Gyr 
with a mass of 3.8 x lO^M© in Salpeter IMF jTanaka et all 



Because most of galaxies with spectroscopic redshift 
are line emitters, a follow-up spectroscopy are needed to 
confirm whether other quiescent galaxies are truly cluster 
members. 

We thus constructed a combined sample of the photo- z 
selected and mass-limited galaxies at 1.56 < z p hot < 1.68, 
composed of 132 [O n] emitters and 259 quiescent galaxies. 
Table 1 gives a summary of our galaxy samples used in this 
paper. 



4 RESULT 

4.1 Spatial distribution 

With the panoramic imaging of the z = 1.62 cluster with 
Suprime-Cam, we have revealed a gigantic structure sur- 
rounding the cluster traced by [O n] emitters for the first 



time, as shown in Figure[7| The cluster appears to be embed- 
ded in a huge filament extending from North to East/South. 
In particular, there are two dense regions to the East/South 
and to the South of the cluster, respectively. Also a filament 
of relatively dense region extends to the North. They all ap- 
pear to be associated to the cluster core and constitute a 
huge structure of about 20 Mpc in comoving scale. 

Many star-forming [O n] emitters are located within the 
projected radius of 1 Mpc in comoving scale from the clus- 
ter centre. This is in stark contrast to the nearby Universe 
where galaxies in dense environment tend to be in-active 
red galaxies and star forming galaxies are preferentially lo- 
cated in lower density environments. We evaluated the over- 
denisity of the [O n] emitters in the cluster core (r<l Mpc). 
The average number density of the [O n] emitters is 0.24 
Mpc" 2 (352/1472) for the entire region of the survey (~1472 
Mpc 2 ). This value is roughly comparable to the number den- 
sity (0.35 Mpc -2 ) expected from th e [O n] luminos ity func- 
tion at z — 1.47 in the general field l|Lv et al.ll2007T ) because 
the criterion of the equivalent width in our selection is more 
stringent than the one used in the previous studies. On the 
other hand, since 13 [O n] emitters exist in the cluster core 
(~3.142 Mpc 2 ), the number density is 4.1 Mpc -2 which is 
about 17 times larger than that in the entire region. The 
quiescent galaxies show the overdensity by factor of 9. This 
clearly indicates that the star formation activity has not 
been ceased yet even in the core of the z = 1.62 cluster, and 
rather the integrated star formation activity in the unit vol- 
ume is actually higher in the cluster core than in the lower 
density regions. 

Note that our narrow-band filter is probing ~94 Mpc 
(comving) in the redshift direction, and the 2-D overden- 
sity may not be free from the projection effect. In order 
to robustly measure the overdensity, the spectroscopic red- 
shifts would be required. In fact, it is found from the FMOS 
spectroscopy that the [O n] emitters in the East and the 
South clumps are located at somewhat different redshifts 
from the cluster centre (Figure 0). This difference corre- 
sponds to 1000-2000 km s -1 in the line of sight velocity 
or to 15-30 Mpc in the comoving distance. It is not clear at 
this stage whether these clumps are gravitationally bound, 
physically associated systems to the cluster. More intensive 
spectroscopic follow-up observations are needed to address 
this issue. 



4.2 Density dependence 

It is widely recognized that the star formation activity of a 
galaxy is strongly related to the local density of galaxies sur- 
rounding it. In the local universe, the star for ming activity 
decreases with increasin g local galaxy density (|G6mez et al.l 
120031 ; iPatel et alj|2009h. This relat i on holds out to at least 
z ~ 0.8 (e.g.. iKodama et all l200ll : iTanaka et alj|2005h . At 
z ~ 1, it is suggested that SFR-density rel ation may, in part , 
be inverted relative to the local relation (|Elbaz et al. I l2007l ; 
ICooper et alj|2008r ). However, these studies do not include 
the very dense environment such as ric h cluster co r es. In 
the C1G J2018-0510 cluster at z = 1.62, iTran et all (|2010l ) 
reported that the relative fraction of star-forming galaxies 
increases with increasing local density based on the SFRs 
derived from the 24 fim fluxes and the SED-fitting. Note 
that they used sample based on only photometric redshift. 
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AR.A. [arcmin] AR.A. [arcmin] 

Figure 7. (Left) The 2-D distributions of 352 [O II] emitters. Blue and red filled circles show star-forming and red [O o] emitters, 
respectively. Squares indicate the spectroscopically confirmed objects with 1.590 < z spC c < 1.620 (black), 1.620< z spC c < 1.630 (green) 
and 1.630 < z Bpcc < 1.644 (red). Contours denote the local number density of [O II] emitters. Outside of dash lines is the field region 
defined in section 4.1. A gray filled circle is a masked region near a bright star. (Right) Close-up view of the cluster core (r< 2'). Dashed 
line circle shows the radius of r=l Mpc. Note that our NB973 filter can prove about ±47 Mpc (comving) in the redshift direction. 



Table 1. The numbers of our photo-z selected samples at 1.56 < 
z p hot < 1-68. 



Sample 


Cluster 


Field 


Total 


star-forming [O II] emitters 


74 


50 


124 


red [O II] emitters 


i 


4 


8 


quiescent galaxies 


86 


173 


259 


combined sample 


164 


227 


391 



On the other hand, iQuadri et all (|2012l l claimed that the 
star formation-density relation holds out to at least z ~ 1.8 
although they similarly selected sample on the basis of the 
photometric redshift. In the case of photo-z selected sam- 
ple, any results must be interpreted with caution due to a 
large contamination. Therefore, the reversal of SFR-density 
relation at z > 1 is very controversial at this stage. 

Our panoramic narrow-band survey of the [O n] emit- 
ters enables us to investigate the environmental depen- 
dence over a much larger area hence covering a much wider 
range in environment. Moreover, due to the great advan- 
tage of narrow-band imaging survey combined with photo- 
metric redshifts, our [O n] emitters are more robust mem- 
bers of star-forming galaxies associated to the cluster even 
without spectroscopy because we just need to separate out 
[O n] emitters at z ~ 1.6 among a few other possibilities 
of different redshifts (corresponding to other lines) based on 
photometric redshifts. In order to discuss the environmental 
effect, we define the local density within the circled area of 
a radius to the 5th nearest neighbour object (Sgth), by us- 
ing the combined sample of z ~ 1.6 galaxies (section 3.3). 
And we calculated the fraction of each sample among the 

combined sample, i.e., / = ATeach sample /Acombimed sample at 



0.8 - 




0.2 - 



0.0 L-d 1 — i 1 I 
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Figure 8. Relative fractions of [O n] emitters (blue circles) and 
quiescent galaxies (red circles) plotted as a function of local den- 
sity. The local density is calculated from the combined samples 
of [O II] emitters and quiescent galaxies. Dashed line shows the 
number density in the field region (~ 1.5 Mpc~ 2 ). In our analysis, 
we did not use the area below this density. Vertical and horizontal 
error bars indicate the Poissonian errors and the sizes of each bin, 
respectively. 



each bin of the local density. However, due to the photomet- 
ric redshift selection, contaminations would dominate the 
number density in the low density regions. Therefore, we do 
not use sample in the field region and concentrated on the 
region within the dashed rectangle in Figure [7] and trust 
only the range of local number density that is more than 
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Figure 9. SFRs (top), stellar masses (middle) and specific SFRs 
(bottom) for individual [O II] emitters, plotted against the local 
density. The local density is calculated from the combined samples 
of [O o] emitters and quiescent galaxies. Red circles show the 
median value in each density bin. 



the field density (Esth, combined > 1-5 Mpc -2 ). The quies- 
cent galaxies sample inevitably contains some foreground or 
background objects because the photometric redshift range 
that we adopted (Az p hot = 0.12) is larger than that of the 
[O n] emitters (Aznb = 0.054). We tried to correct for this 
effect by multiplying Aznb/ Az p hot to the number of qui- 
escent galaxies, assuming conservatively that their redshift 
distribution is uniform. As shown in Figure[8l the fraction of 
star-forming galaxies in the combined sample do not shows 
a significant dependence on the local density at z ~ 1.6. 
It is consistent with a constant value (~60%) across differ- 
ent environments within errors. Star-forming activity in the 
cluster core is very high in this cluster, and the well estab- 
lished star formation-density relation in the local Universe 
no longer exists. Although a lot of massive, quiescent galax- 
ies do exist in the high density environment such as the 
cluster core, the star formation activity has not been ceased 
yet. Rather the integrated star formation rate per unit vol- 
ume is actually much higher in the cluster core due to its 
high over-density of the emitters. 

In order to quantitatively compare the properties of 
the individual [O n] emitters as a function of environ- 
ment, we estimate stellar mass, SFR and specific SFR 
(SSFR=SFR/M st ar) for the [O n] emitters. [O n] emission 
line is widely used as a good SFR indicator but strongly 
depends on metallicity and dust-extinction. [O n] -SFR re- 
lation ar e recently being calibrated by using low er redshift 
samples (|Gilbank et al.l l2010l ; ISobral et al.ll2012i). We use d 
the [O n] -SFR calibration given by ISobral et all (|2012i ). 
which is described as a function of stellar mass. The results 
are shown in Figure [5] as a function of local density. We do 
not find any significant environmental dependence in all of 
the individual quantities of the [O n] emitters. In fact, the 
median values of a ll the three quantitie s are almost constant 
with local density. iTadaki et al.l |201ll ) claimed that the av- 
eraged star formation activity is rapidly declined from z=2.5 



to z=0.8 in high density regions such as clusters while it is 
only gradually declined in the low density regions. At some 
point at high redshifts, we may expect that the averaged star 
formation rate of galaxies in cluster cores should exceed that 
in the general field probably because galaxy formation pro- 
cesses (such as gas cooling and mergers) are accelerated in 
dense environments. The fact that we observe just a com- 
parable level of star formation activity in the [O n] emitters 
irrespective of environment may suggest that, at z = 1.62, 
the [O n] emitters in the high density regions are just in the 
transition phase from a star-bursting mode to a quiescent 
mode as a result of environmental effect. 



5 SUMMARY 

We have conducted a narrow-band survey of [O n] emit- 
ters in and around the C1G J0218. 3-0510 cluster at z = 
1.6, using Suprime-Cam on Subaru. The observation with 
z-r filter was newly carried out to measure the contin- 
uum at the same wavelength as the narrow-band filter 
(NB973). We combined these data with multi broad-band 
data (B, V, R, i' , z' , J, H, K, 3.6 /mi, 4. 5 fim) to create the 
NB973-detected and K-detected photometric catalogues. 
The photometric redshifts were determined with 11-bands 
including z K -band. By comparing the spectroscopic mem- 
bers and our [O n] emitter sample, we estimated that the 
accuracy of our photometric redshifts is <j 2 = 0.05 at z = 1.6. 
Our main results with these samples are summarized as fol- 
lows. 

(i) On the basis of narrow-band excesses and photometric 
redshifts, our survey provides a sample of 352 [O n] emit- 
ters over a 830 arcmin 2 area. Our very recent FMOS 
near-infrared spectroscopic observations have confirmed 31 
[O n] emitters at z ~ 1.6 by the presence of Ha or [O ill] lines 
at the expected wavelengths. The [O n] emitters constitute a 
large scale structure at z = 1.62 in which the C1G J0218.3- 
0510 cluster is embedded. Also, we find that many star- 
forming [O n] emitters are located even in the cluster core 
(r<l Mpc in the comoving scale) and in the surrounding 
clumps, and show a high overdensity by a factor of 17 com- 
pared to the entire region. This suggests that the integrated 
star formation activity per unit volume is activated in such 
regions. 

(ii) The galaxies with 1.56 < z pa ot < 1-68 show a clear 
bimodal distribution in the zJK diagram, namely, quiescent 
and star-forming galaxies. We selected total of 259 quiescent 
galaxies on the basis of the population synthesis model and 
the extinction law by dust. 

(iii) We calculated the local number density of galaxies 
in the photo- z selected and mass-limited sample to examine 
the environmental dependence of the star-formation activ- 
ity. We obtain a large fraction of [O n] emitters even in the 
cluster core, showing that the star forming activity in the 
cluster core is elevated substantially compared to the local 
clusters where there are little star-forming galaxies in their 
cores. There is no longer a environmental dependence in the 
relative fraction of [O n] emitters in the combined sample, 
and the well known SFR-density relation in the present- 
day Universe no longer exists within errors. Furthermore, 
the properties of the individual [O n] emitters, such as star 



10 K. Tadaki et al. 



formation rates, stellar masses and specific star formation 
rates, do not depend on the local density, either. These re- 
sults suggest that the [O n] emitters in the high density 
regions are just in the transition phase from a star-bursting 
mode to a quiescent mode due to some environmental ef- 
fects and the star formation rates in these systems may be 
rapidly declining. 

Note that the results presented in this paper for the en- 
vironmental effect of galaxies in and around a 2=1.62 clus- 
ter may not necessarily apply to all the clusters at the same 
epoch of the Universe, and the results may well depend on 
the degree of matureness of clusters even at the same red- 
shift. For this reason, we require a systematic studies of dis- 
tant clusters and general fields to construct a more general 
picture, and the MAHALO-Subaru project (Kodama et al. 
in prep) will provide us more comprehensive views of galaxy 
evolution at its most active phase in the Universe. 
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